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Selection of dune shapes and velocities. 

Part 1: Dynamics of sand, wind and barchans. 
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Abstract. Almost fifty years of investigations of barchan dunes morphology and dynamics is reviewed, 
with emphasis on the physical understanding of these objects. The characteristics measured on the field 
(shape, size, velocity) and the physical problems they rise are presented. Then, we review the dynamical 
mechanisms explaining the formation and the propagation of dunes. In particular a complete and original 
approach of the sand transport over a flat sand bed is proposed and discussed. We conclude on open 
problems by outlining future research directions. 

PACS. 45.70.-n Granular systems - 47.54.+r Pattern selection; pattern formation 



1 Introduction 

After pioneering works by Bagnolds Q, the investigation 
of dune morphology and dynamics has been the exclu- 
siveness of geologists and geographists. For nearly four 
decades, they have been reporting field observations about 
the condition under which the different kinds of dunes ap- 
pear as well as measuring their velocity their 
shape the patterns they form, 

the size distribution of sand grains, etc. More recently 
dunes have attracted the interest of physicists ]p^,[T3|,p^, 
[T^,|l6,r7,l^,[l9| motivated by the poetry of deserts or by 
the physics they conceal, by the dryness of saharan coun- 
tries or by the observation of dunes on Mars. 

The aim of this paper is to give an overview of this 
domain. It is conceived as a pedagogical review of field 



observations and of the dynamical mechanisms impor- 
tant for dune morphogenesis, followed by a conclusion 
discussing some of the problems remaining open. It is fol- 
lowed by a second part devoted to the derivation of a 
simple model predicting the shape and velocity of two- 
dimensional dunes. 

In this article, we will mainly focus on the barchan 
which is the simplest, and consequently the most studied 
form of sand dunes. A barchan is a crcscentic dune as that 
shown on figure ^, propagating downwind on a firm soil 
[0,|2^]. When the direction of the wind is almost constant, 
these dunes can maintain a nearly invariant shape and size 
(3— 10 m) for very long times (1—30 years) The basic 
dynamical mechanism explaining the dune propagation is 




Fig. 1. The barchan is a dune with a characteristic crescent 
shape, possibly isolated. A small dune (~ 3 m high) as that 
shown (Mauritania) , propagates downwind at one hundred me- 
ters per year, typically. 



Fig. 2. Sketch of a barchan dune. In first approximation, the 
dune morphology can be described by four morphologic pa- 
rameters: the length L, the width W, the height H and the 
horn length Lhorns- 
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simple (see figure ||) : the back of the dune is eroded by the 
wind; the sand transported in saltation is deposed at the 
brink, and is redistributed on the slip- face by avalanches. 
In this picture, a dune keeps the same amount of sand. In 
reality, the global mass evolution results from the balance 
between the sand supply at the back and the leak at the 
horns. 

Starting from this overview of barchan properties, we 
give in section ^| a short review of the field observations, 
concerning in particular the morphology and velocity scal- 
ing laws. The dynamical mechanisms involved in erosion 
and sand transport are presented in section |3| and the cor- 
responding dimensionless parameters as well as the scal- 
ing laws are discussed. Finally, in section ^, we sort the 
facts both validated by experimental measurements and 
explained by consistent theories from the problems re- 
maining partly or totally open up to now. 



2 Field observations 



L (m) 
80- 



60 — 



40 



20- 



o- 




»Finla-l lOW (44 h>m :h>ms ) 
♦ fjjiig $ Sharp 1 9IS4 (47 harrhfmO 
Al [aSUSriCMti 1967 (42 lawliaii*) 
AlijKkiirdtli 1987 (6 basilar*} 







H (nO 



Fig. 3. Relationship between the dune length L and height 
H determined from field measurements averaged by ranges of 
heights. The solid line is the best linear fit to the points cor- 
responding to barchans from the Arequipa region in Southern 
Peru (Finkel M and Hastenrath 
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Fig. 4. Relationship between barchan width W and height 
H determined from field measurements averaged by ranges of 
heights. The solid line is the best linear fit to the points cor- 
responding to barchans from the Arequipa region in Southern 
Peru (Finkel § and Hastenrath [|,|). 



2.1 Morphologic relationships and scale invariance 

The crescent-like shape of the barchan is well known (see 
figures l] and ^|) and can be characterized by a few param- 
eters, the length L along the central axis (which is also the 
wind direction), the height H, the width W and the horn 
length Lhoms (see figure ||). Note that in practice, the two 
horns have always different lengths due to the fluctuations 
of wind direction. They are thus measured separately and 
averaged to give L horns . 

The region around La Joya in southern Peru is the 
most documented barchan field and that for 

which measurements are the most coherent. Finkel [[| and 
later Hastenrath ||^| have chosen a squared region, a few 
kilometres width, and have investigated systematically the 
morphologic parameters of the barchans in this perimeter. 
Their data exhibit a large statistical dispersion due to the 
variations of the control parameters on the field. Then, it 
is not surprising to find a dependence of the dune shape on 
the local conditions like the wind regime, the sand supply, 
the presence or not of other dunes in the vicinity, the na- 
ture of the soil, or wether the studied dunes have achieved 
a permanent state or are in a transient state. . . 
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Despite this strong dispersion of data, clear linear re- 
lationships between the height H, the length L and the 
width W (see figure ||) were found by Finkel (|] and Has- 
tenrath B, but with different coefficients (W = Wo+pwH 
and L = Lo + plH). To seek for relationships valid sta- 
tistically, we follow the procedure of Finkel Q and aver- 
age all the measurements by ranges of heights (see fig- 
ures H and More precisely, we have considered the 
barchans between 1 to and 2 m high and averaged their 
length, width and height; this gives one point on figures [| 
and ^. The same computation was done for the barchans 
between 1.5 m and 2.5 to, between 2 to and 3 to, etc. 
We can reasonably expect the statistical relations to be 
that which would have been obtained in controlled, repro- 
ducible conditions. 

It can be observed that linear relationships between 
the length L, the width W and the height H are nonethe- 
less recovered but are actually the same for the three sets 
of measurements (Finkel J3| and Hastenrath made 
in Southern Peru. The best linear fits give pw = 8.6 
and pl — 5.5 corresponding to the mean proportions of 
barchans near La Joya. The most important point is that 
L, H and W are not proportional. On the contrary, the co- 
efficients L = 10.8 to and Wo = 8.8 to are comparable to 
the size of the smallest dunes. This means that small and 
large dunes are not homothetic: dunes of different heights 
have different shapes. In other words, barchans are not 
scale invariant objects. As a consequence, there should ex- 
ist (at least) a typical lengthscale (related to Lo and Wo) 
in the mechanisms leading to the dune propagation. 

There are very few places in the world similar to the 
Arequipa region (southern Peru) . They are usually touched 
by trade winds driven by oceanic anticyclones. Because 
these anticyclones are very stable, the winds are constant 
and often very strong. It is the case with the Peruvian 
and Chilean coasts, with coastal Namibia, with the At- 
lantic coast of the Sahara from Morocco to Senegal and 
with the northern shores of Western Australia. Morpho- 
logic measurements were conducted in some of the places 
where there is a constantly prevailing wind direction and 
a smooth ground surface: in Mauritania by Coursin 
in California by Long and Sharp Is], and more recently in 
Namibia by Hesp and Hastings pi and in southern Mo- 
rocco (former Spanish Sahara) by Sauermann et al. ]l2|] . 
In the later case, the authors have not measured the whole 
barchans in a given area but selected some which where 
isolated and symmetric. Unfortunately, in these works, the 
number of barchans is too small to determine precise mor- 
phologic relationships. But they are sufficient to indicate 
that the morphology depends on the dune field (see fig- 
ures H and ||). Briefly, we know almost nothing on the pa- 
rameters important for dunes morphology. It can depend 
on the grain size, the density of dunes, the sand supply, 
the wind strength, its changes of direction, etc. 

For instance, it has been observed by Allen [^3| that 
barchans have a smaller width W and more developed 
horns Lhoms under strong winds than light ones. But 
Hastenrath have reported measurements which indicate 
a weak dependence of the shape on the wind strength. He 



returned to southern Peru |Q] after a light wind decade 
and observed that all the dunes had strongly decreased in 
size. He measured a second time the morphologic param- 
eters of 6 barchans (dotted triangle on figures || and ||) : 
the points are almost on the same line than previously. 

Another example is the coincidence (or not) of the 
brink and the crest (see figure || for a definition of these 
words). Hastenrath || and Sauermann et al. [l2]| have ob- 
served that small dunes present a broad domed convexity 
around the crest, clearly separated from the brink as in 
figure |2|, while large dunes have the crest straight to the 
brink as in figure 0. This is represented schematically on 
figure |Io|. However, it had been reported before [^(j) that 
there exist barchans of the same height in the same dune 
field presenting alternatively the separation or the coinci- 
dence of the brink and the crest. The difference lies per- 
haps in the selection of barchans to be studied (isolated 
or not, symmetrical or not, etc). 

The last morphologic relation between the mean length 
of the horns Lhoms and the barchan height H (figure ||) 
is a nearly perfect relation of proportionality: Lhoms — 
9.1 H. Even more striking, the relation is approximately 
the same for all the dune fields measured (in particular La 
Joya, Peru an d Imperial Valley, California J|], see 

figure H). This suggests that the scaling law of the horns 
length could be simpler and more robust than that of the 
back dimensions. 
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Fig. 5. Relationships between the horn length Lhoms and 
the width W determined from field measurements averaged 
by ranges of heights. The solid line is the best linear fit to the 
points corresponding to barchans from the Arequipa region in 
Southern Peru (Finkel Q and Hastenrath |^[§). 
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Fig. 7. Time evolution of a sandpile blown out by a controlled wind without (a) and with (b) a sand supply. In both cases 
the sandpile is eroded and disappears after a few minutes. Without sand supply, the erosion is localised on the sides so that a 
longitudinal brink is created. With a sand supply, the pile takes the form of a crescent. In particular, the back remains smooth 
and two horns grow by deposition of grains in reptation. In case b), both top an side views are shown. In particular, the initial 
height of the pile (~ 20 cm) can be seen. 
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Fig. 6. To study the time evolution of a sandpile blown by the 
wind, the Cemagref wind tunnel (6 m long, 1 m large and 1 m 
high) was used. Pictures were taken from above using a video 
camera. A mirror was placed at 45° to get on the same picture 
a side view and a top view of the pile. The later was enlightened 
by a lamp and by an horizontal laser sheet adjustable in height 
which reveals the topography. A tunable sand supply has been 
add at the top of the tunnel. The grains are PVC beads of 
size 100 (jm. The velocity (around 6 m/s at 2 cm above the 
soil covered by velvet) is chosen slightly above the threshold of 
motion of the grains. 



2.2 Minimum size 

It is striking to note the absence, on figures pi |] and 0, 
of measured dunes smaller than H = 1 m, W = 19 m 
and L — 17.5 to. This cut-off is clearly visible on the dune 
size histograms measured by Hastenrath ||. Moreover, if 
barchans can be very thin sand patches in the region where 
they form, they are never lower than, say, 1 to in mature 
dunes regions. What appends to a small barchan ? Re- 
lated to this question, there have been several attempts 
|jl|,|o[|2^] to generate an artificial dune from a small coni- 
cal sandpile (typically 10 cm to 1 to high). 



Here we report an experiment made in the Cemagref 
wind tunnel, described on figure ^. A conical sandpile 
20 cm high is built from a funnel. It is then eroded by the 
air flow with (case b) or without (case a) a sand supply 
at the beginning of the wind tunnel (figure 0). Whatever 
the conditions (even with a sand supply), the pile looses 
mass inexorably, and disappears in a few minutes. This is 
also what was found in field experiments JH|20|,^l[. 

Figure shows the evolution of the sandpile in the 
two cases, for a wind velocity chosen slightly above the 
threshold of motion of the grains. In case a), without a 
sand supply, the erosion is localised on the sides and the 
formation of a longitudinal brink is observed. In the last 
steps the side faces become so steep that avalanches occur. 
In case b) , the input sand flux is tuned to be slightly below 
the saturated flux. Both the back and the sides are eroded 
so that the pile takes the form of a barchan. In particular, 
we observe the progressive formation of two horns due to 
the lateral deposition of grains in reptation. 

The observation that a sandpile 20 cm high disappears 
whatever the conditions suggests that barchan dunes have 
a minimal size. This was first noticed by Bagnolds who in- 
terpreted the cut-off scale (the minimum dune size) as the 
saturation lengt h wh ich will be defined and discussed in 
details in section 3.3. Basically, it is the length over which 
the sand carrying increases when the wind passes from 
a firm soil to a sand patch. This interpretation has not 
been confirmed so far and deserves further investigations. 
Still the existence of a minimal size rises crucial problems. 
First, it means that no small size dune can be obtained in 
the air. Second, it asks a fundamental question : if a small 
barchan (or a conical small sandpile) is unstable and dis- 
appears after a short time, how then can barchans form? 



2.3 Barchans velocity 

In marryplaces, barchans velocity has been measured |], 
|,||,§|Jl0). Typically, a small dune, 3 to high, propa- 
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gates at a speed ranging from 15 m/year to 60 m/year 
while this velocity is between 4 m/year and 15 m/year for 
a large dune, say 15 m high. Even though measurement 
points are very dispersed |2Q] , there is no doubt that small 
dunes move faster than large ones. 
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Fig. 8. Dune velocity, averaged by ranges of height, as a func- 
tion of height. Solid lines correspond to the best fit of the data 
by a velocity of the form c = Q/(Hq + H). The fit is shown 
for the measurements in the Arequipa region (Finkel |J and 
Hastenrath [0,0]) 



As for morphologic relationships, we have grouped the 
measurement points and averaged the velocity over several 
barchans. The resulting curves are shown on figure |^. It 
can be seen that the dune velocity is very different from 
one place to another. In the imperial valley, California, 
dune displacements were measured by Long & Sharp H 
from 1941 to 1956 (black diamonds), and from 1956 to 
1963 (white diamonds). In Peru, they were measured by 
Hastenrath ||||] from 1955 to 1958 (black triangle), from 
1958 to 1964 (white triangle), and from 1964 to 1983 (dot- 
ted triangle). From these two data sets, it can be seen 
that barchan velocity c strongly depends on time. This 
is probably related to the fact that it obviously depends 
on fluctuating parameters like the wind speed or the sand 
supply. 

However, in each dune field and especially in Southern 
Peru (figure ^|) , the relationship between velocity and size 
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Fig. 9. Inverse dune velocity 1/c, averaged by ranges of height, 
as a function of height H . Only the measurements made in the 
region of La Joya (Southern Peru) are shown. 



can be reasonably described by: 

Q 



H Q + H 



(1) 



This relation is represented schematically on figure |T(j. Q, 
which is homogeneous to a volumic sand flux, is found 
to be approximately equal to 85 m 2 /year between 1955 
and 1958 and to 425 m 2 /year between 1958 and 1964, in 
southern Peru (figure ||). Hq is a cut-off height ranging 
from 1.8 m (Finkel 1955-58) to 10.9 m (Hastenrath 1958- 
64). 

^From the propagation speed, we can construct the 
dune turnover time T turn over as the time taken by the 
dune to travel over its own length: 



L _ L(H Q + H) 
c~ Q 



(2) 



It is also the typical period of the cycle of motion of a 
grain of the dune: erosion from the back, deposition at 
the top of the slip face, flowing in an avalanche, at rest 
below the dune, reappearance at the back of the dune, 
etc. It gives also the typical time needed by the dune to 
readjust its whole shape to changes of external conditions 
(wind, sand supply. . . ). Typically, a small dune, 3 m high, 
has a turnover time between 5 months and 2 years while 
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Fig. 10. Visual representation of the relation between the dune velocity c and its height H: three barchans are represented at 
initial time and after some time t. It can be seen that small dunes go faster than large ones. From the scalings of the speed c 
and the width W (see text), we can infer that the area W ct swept by the barchan (in light grey) is almost independent of the 
dune size. It can be seen on the large barchan that the back of the dune is eroded while the slip face and the horns are regions 
of sand deposition. 



it ranges from 6 years to 25 years for a large dune, 15 m 
high. According to Oulehri ^4j, this large difference of 
turnover times between small and large dunes was already 
known in ancient times by saharan people. Barchans were 
used as cereal lofts or to protect goods from pillaging. A 
small or a large dune was chosen, function of the time 
after which they wanted to recover the bundle. 

3 Mechanisms, dimensionless parameters and 
scaling laws 

The previous field observations and measurements rise 
several questions. What are the basic dynamical mech- 
anisms acting in the initiation and the propagation of 
barchans? What determines the equilibrium shape of bar- 
chans? Can we predict the morphologic relationships and 
the speed of dunes? Why is barchan shape not scale in- 
variant? 

The first part of this program i.e. the basic dynamical 
mechanisms, have been investigated in Bagnolds pioneer- 
ing work ||]. We give here a complete overview of these 
mechanisms and we reformulate the corresponding scaling 
laws. 



3.1.1 Turbulent boundary layer 

The dune dynamics is controlled by the sand transport 
which is itself driven by the wind. For a fully developed 
turbulent wind over a flat surface composed by grains of 
typical size d, the wind velocity u usually increases log- 
arithmically with height z (see figure [ll]). This can be 
simply understood in the context of turbulent boundary 
layers theory. The standard turbulent closure relates the 
air shear stress r to the velocity gradient d z u: 

where p ,> is the volumic mass of air (see table and n ~ 
0.4 is the Von Karman constant. For a steady and uniform 
boundary layer, the shear stress r is constant and equal 
to To, the shear force per unit area on the bed. Besides, 
it turns out that the velocity vanishes at a distance rd 
from the sand bed. The rescaled bed roughness r is found 
to be of the order of 1/30 Jl],|o|. This gives a logarithmic 
velocity profile, as observed on the field: 



3.1 Saltation and reptation 



where the shear velocity u* is by definition u* = \Jtq/ p air . 
We will see in the following how the sand transport mod- 
ifies this velocity profile. 
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3.1.2 Motion of one grain in the wind 

A wind of sufficient strength can dislodge and entrain sand 
grains. Once they have taken off the sand bed, they are 
progressively accelerated by the wind, as they go up. The 
trajectories are observed to be asymmetrical, as shown on 
figure |ll|. The grains are submitted to the gravity and to 
the fluid drag force which leads to the standard following 
equation: 



rfv 
~dt 



= g + X 



Psand 



(5) 



where v is the grain velocity, p sa nd the volumic mass and 
u the local wind velocity - which depends on z. The drag 
coefficient \ depends on the shape of the grain but also 
on the Reynolds numbers Re = d |u — v\/v {y stands for 
the air viscosity). At large Re, in the turbulent regime, \ 
tends towards a constant if the grain is sufficiently rough 
p7[. In the viscous - Stokes - regime, \ decreases as 1/Re. 




Fig. 11. In a fully turbulent boundary layer, the velocity pro- 
file is logarithmic. A grain dislodged from the bed is accelerated 
by the wind as it goes higher and higher and decelerated in the 
descent so that its typical trajectory is asymmetric. Saltation 
hop lengths l sa i are found to be about 12 to 15 times the height 
of bounce Z sa i ■ 



In the vertical direction, the motion is dominated by 
gravity. It immediately follows, as for a free parabolic 
flight, that the flight time T and the height of bounce 
h depend on the launch speed w as: 



1 oc — ; 
9 



Z oc — . 
9 



(0) 



Note that, due to the quadratic drag force, these expres- 
sions must be corrected at strong winds |^5| : the height of 
bounce Z get actually smaller than that of the free flight. 
These scalings are valid whatever the launch velocity w 
is. In the sequel we will distinguish between the different 
kinds of trajectories. For example, the indexes 'sal' and 
'rep' will be used for saltation and reptation grains. If no 
label is specified, it means that the argument applies in 
general. 

How does the flight time T compare with the timescale 
Tdrag after which the grain has been accelerated by the 



drag force up to the wind velocity? Using the equation of 
motion (||), we get: 



T, 



X Pa 



(7) 



For typical values of the parameters (table [j]), the saltation 
flight time T sa i and the drag timescale Tdrag turn out to be 
of the same order of magnitude. This means that grains - 
larger than 100/xm - are not much sensitive to the turbu- 
lent fluctuations of the wind but also that the horizontal 
velocity of the grain is not far from the mean wind velocity 
along the trajectory. The grain spends most of the flight 
time around the trajectory maximum (z ~ Z). As a con- 
sequence, the grains mean horizontal velocity should be 
approximately equal to the wind velocity u at this height. 
The horizontal displacement - the saltation length - is of 
the order of the flight time T times the mean horizontal 
velocity: 

ran 

l(X ~ • ( 8 ) 



28 



and indirect [£6,27 measurements 
29(1 are in agreement with this 



All the direct 0. 
as well as the models 
simple approach of grains trajectories. On the other hand, 
the scalings of Z and / with the wind shear velocity it* 
and the grain diameter d are still controversial up to now. 
Mainly, the problems are the modification of the wind - 
and thus u - by the saltating grains and the mechanism 
by which the launch velocity w is selected. As from now, 
it can be argued that the fastest grains, which are said 
to be in saltation, bounce so high that the wind u at the 
trajectory maximum Z sa i is almost undisturbed by the 
rare grains to pass there. Their velocity, in particular the 
vertical component w after a collision, should thus scale as 
u„ . Taking w = m» , as suggested by Owen Q , the height 
of bounce reads: 

it 2 

Z snl ~ (9) 



The horizontal velocity u (at z = Z. 
with a non dimensional prefactor (u 
logarithmic velocity profile: 



sai) scales with it* but 
: £ u*) reflecting the 



« l \r 9 d 



For those high energy grains, the saltation length, 



9 



(10) 



(11) 



is £ times larger than the saltation height Z sa i (figure "Tj) . 
Experimentally Jl|,^5), the ratio l sa i/Z sa i is between 12 
and 15, which is the order of magnitude found here for 
£ (table 0). The grains bouncing the highest should thus 
have trajectories almost independent of their size d and 
scaling on ti^/j for the lengths and onii, for the velocities. 
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Table 1. Definition and typical values of the main quantities 
discussed in the text. 



Dune length 


L 




Dune height 


H 




Dune width 


w 




Length of dune horns 


Lhorns 




Dune velocity 


C 




Dune profile 


h(x,y,t) 




Sand flux 






Rescaled roughness 


r 


1/30 


Von Karman constant 


K 


0.4 


Shear stress 


T 




Air-borne shear stress 


Tair 




Grain-borne shear stress 


7~sand 




Quartz density 


Psand 


2650 kg.m' 3 


Grain diameter 


d 


200 urn 


Shear velocity 


It* 


0.5 m.s' 1 


Logarithmic velocity prefactor 


s 


20 


Wind velocity 


u 


10 m.s" 1 


Saltation flight time 


T S al 


50 ms 


Saltation hop height 


Z S al 


2.5 cm 


Saltation hop length 


^sal 


50 cm 


Reptation hop length 


I rep 


8 mm 


Static friction coefficient 


Ms 


tan (33°) ~ 0.65 


Dynamic friction coefficient 


fid 


tan(31°) ~ 0.60 


Air density 


p fluid 


1.2 kg.m' 3 


Air viscosity 


V 


1.5 10" 5 m 2 .s~ 


Grain Reynolds number 


Re 


du»/v ~ 120 






du/v ~ 2500 


Impact threshold velocity 




15 cm.s -1 


Fluid threshold velocity 


Uflu 


20 cm.s' 1 


Saturated sand flux 


Qsat 


180 m 2 /year 


Drag time 


Td ra g 


50 ms 


Drag length 


^drag 


9 m 


Reptation flux 


Qrep 




Vertical sand flux 






Number of splashed grains 


Neje 





Poo ~ 0.95 is the rebound probability for velocities much 
larger than a^/gd. a is a non-dimensional number equal to 
10 in Anderson and Haff two-dimensional simulations [|TJ 
EH- 

Experiments ]25|,[33| and numerical simulations [ ^ . [3^ ] 
show that the rebound velocity fraction of the 



impact velocity v^, 



Vreb = 7 V i 



The restitution coef- 



ficient 7 is around 0.5. The rebound angle Q re \, is almost 
independent of the impact velocity (modulus and angle) 
and ranges from 35° to 50°. 




Fig. 12. When a saltating grain collides the sand bed, it re- 
bounds and splashes up other grains. 



The characteristics of the ejecta are also well docu- 
mented [^l[|3^,|3^| . The distribution of the ejection veloc- 
ity - both the modulus v e j e and the angle 9 e j e from the 
horizontal - is found to be almost independent of the im- 
pact velocity Vi mp . The mean ejection speed is related to 
the escape velocity y/gd: 



a-\fgd and e 



70° 



(13) 



3.1.3 Collision of one grain on the sand bed 



A constant fraction of the impact momentum is trans- 
ferred to the ejecta so that, on the average, the number of 
ejecta increases linearly with the impact speed: 



When a saltating grain collides the sand bed, it rebounds 
but can also eject other grains (figure [jj]). In general, 
the rebounding grain differs sufficiently from the ejected 
grains to be recognised. However, there are low probability 
configurations for which the incident grain delivers its mo- 
mentum to few nearby grains but remains trapped, even 
for large impact velocities. The non dimensional parame- 
ter which controls the rebound probability p re b is the ratio 
of the impact velocity Vi mp to the velocity necessary to es- 
cape from the potential trapping at the sand bed surface 
[p0|| , namely \fgd. In particular p re b vanishes when Vi mp 
becomes smaller than this escape velocity. The rebound 
probability resulting from the numerical simulation of An- 
derson and Haff |3llB2] can be expressed as: 



Preb — Pc 



exp 



J imp 



a^/gd 



(12) 



N eje = 



y imp 



ay/gd 



1 



if Vimp ^> CL 



gd 



(14) 



Again, the same critical velocity a^fgd appears. It is at 
the same time the critical velocity below which no grain 
is on the average ejected, the mean ejection speed and 
the velocity below which the rebound probability strongly 
decreases. 

These low energy grains ejected from the sand bed 
move near the surface of the sand bed: they bounce typ- 
ically at a few hundred times the grain size. Such grains 
are said to be in reptation J34|, in contrast with saltation 
which relates to high energy bouncing grains transported 
by the wind. The reptation is at the origin of the forma- 
tion of ripples which, by the way, propagate and thus also 
takes part in the sand transport. 
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3.1.4 Saltation versus reptation 



The slowest grains, the 'reptons', have properties that de- 
pend only on the grain size: the bouncing height scales 
with d and the grain velocity with \fgd. On the other 
hand, the fastest grains, the 'saltans', have properties which 
depend mainly on the shear velocity it*. Between the two, 
there is a continuum of trajectories of different jump lengths 
which all contribute to the sand flux (3^, ^,^9). Experi- 
mentally, this can be observed directly or deduced from 
the fact that the sand transport is a continuous function 
of height. Why then making a distinction between repta- 
tion and saltation? 

The first answer is of course that the two limit trajec- 
tories - slowest and fastest - of the distribution do not 
present the same sc alin g. Let us turn to the experiment 
presented in section |2.2| (figures ^| and [?]) which clearly il- 
lustrates the difference. For any wind velocity above the 
threshold of motion of the grains, the saltation length is 
larger than the size of the pile (typically 1 m compared 
to 40 cm): no saltating grain dislodged from the pile is 
deposed on it. As a consequence, the pile is only eroded. 
Then it is surprising to observe on figure ^ the formation 
of horns, downwind the initial position of the pile. This 
means that they have been formed by grains deposed there 
which can only be grains in reptation. The whole figure [7] 
can be read in this way: the pile is eroded due to saltation 
but the formation of a crescent shape in particular the 
horns are due to reptation. In this situation, the flux of 
grains in saltation (visible by the overall leak of matter) 
and the flux of grains in reptation (visible by the forma- 
tion of horns) are comparable. 

and numeri- 



As indicated by experiments 136,3^,^ 
cal simulations [ p^l{|32]] the sand transport is maximum at 
the ground level and decreases exponentially with height. 
Thus, the contribution of reptons to the overall flux of 
sand should be important. The ratio between saltation 
and reptation fluxes has only been measured directly by 
Bagnolds Q using two different traps: a rectangular trap 
was placed vertically for the saltons and a second one, 
with a thin linear mouth, was buried in the soil to trap 
the grains moving just at the surface. In fact this hole 
could only trap the slowest and largest reptons. Bagnolds 
proposed to designate this motion by surface creep. An- 
derson, S0rensen and Willetts Q gave a more precise def- 
inition of creeping grains as grains which get rearranged 
by saltation impacts and which are not affected directly 
by wind forces. This category is actually useful for binary 
mixtures Q, to describe the motion of the heavy grains 
submitted to a rain of light ones. With these traps, Bag- 
nolds found that the creeping flux was only three to four 
times smaller than the saltation flux. 

Since the grains in reptation are dragged along by 
saltating grains, we can infer that the reptation flux is 
just proportional to the saltation flux. Defining the verti- 
cal flux <f> as the volume of particles leaving the soil per 
unit time and unit area, the reptation flux <j) rep is governed 
by the splash of saltation grains. In first approximation, 




Fig. 13. The flux q is the volume of sand which crosses a unit 
line transverse to the wind per unit time. If the typical path of 
a grain has a hop length I, the incident flux of grains <fr, which 
is the volume of grains colliding a unit area per unit time, is 
equal to q/l. 



this can be written as: 



N e je <t>sal 



(15) 



Because N e j e increases linearly with u^/y/gE, the number 
of reptons leaving the soil per unit time is much larger 
than the number of saltons. For the sand transport, the 
important quantity is the flux q defined as the volume of 
sand crossing a unit line transverse to the wind per unit 
time (figure 13). The relation between the vertical flux <fi 
and the integrated horizontal flux q is determined by the 
grains typical hop length I. For the grains in saltation, this 
gives: 

qsal = 4>sal Lai (16) 

Then the relation ( filf ) between fluxes of saltation and rep- 
tation can be expressed in terms of q: 



Irep N e j e 

Qrep — 1 Qsal 
^sal 



(17) 



Since l re p/hai decreases as gd/u*, this means that the con- 
tribution of reptons to the overall flux of sand decreases 
with the wind velocity. This is consistent with the disap- 
pearance of ripples at large wind. One should thus be very 
careful in identifying the species most important for the 
sand flux: the result is opposite for the horizontal flux q 
and the vertical flux <j>. Finally, the same question can be 
risen: if the two fluxes are proportional one with the other, 
why making a distinction between the two? 

Another important difference appears when the sand 
bed presents a slope, in particular if the steepest slope 
is perpendicular to the wind. Howard Q have observed 
that the normal to ripples are no more parallel to the 
wind on sloping surfaces but are deflected downslope by 
as much as 35°. This clearly shows that the gravity force 
has an influence on the direction of motion of reptons. 
Moreover, Howard |l0| found a good agreement with a 
simple calculation considering that the ejected grains are 
submitted to an effective drag along the wind direction, 
just sufficient to escape from potential trapping. Most of 
their flight, saltons are dragged along the wind direction 
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but their rebound direction could be influenced by the 
surface slope. 

More recently, Hardisty and Whitehouse EIJ] have shown 
that the bedslope along the direction of motion modifies 
the threshold of motion (typically 15% for a 10° slope) but 
also the flux itself. They found that the flux was multi- 
plied by 10 for a 10° downslope and divided by 6 for a 10° 
upslope. These very strikin g re sults were not confirmed 
by Rasmussen and Iversen |3^,[39],[l2| who have found a 
much weaker dependence of the saturated saltation flux 
with the longitudinal slope. 

There still misses experimental studies investigating 
the dependence of the reptation flux with the slope - in 
particular a lateral slope. Still, it can be argued that grav- 
ity is important in the ejection process. The number of 
ejecta should be smaller for a positive slope along the wind 
direction than for a negative slope and reptons should be 
deflected downslope by gravity. Then, the total flux is not 
aligned along the direction of wind but has a component 
along the steepest slope. To the first order this downs- 
lope flux is simply proportional to the gradient V/i of the 
quantity transported (the local height h) . It thus leads to 
diffusive effects which and tend to smooth the dune. 



3.2 Entrainment by the fluid and by impacts 

3.2.1 Fluid threshold 



Table 2. Typical values of the quantities denned and discussed 
in the text. The aeolian sand transport is compared to that in 
water. 




Fig. 14. Wind is accelerated over grains at the surface of the 
sand bed. To lift up a grain, the low pressure created by the 
wind flow has to be larger than gravity. There is thus a thresh- 
old velocity below which the wind can not erode the sand bed. 



wind: 



Fluid 


Air 


Water 


P fluid 


1.2 kg.m' 3 


1000 kg.m,- 3 


V 


1.5 10~ 5 m 2 .s" 


1 lO" 6 m 2 .*" 1 


Re, = du*/v 


120 


8 


Re = du/v 


2500 


170 


Uimp 


15 cm.s -1 


15 cm.s- 1 


Uflu 


20 cm.s" 1 


0.7 cm.s -1 


q S at 


180 m 2 /year 


17 m 2 /hour 


Tdrag 


50 ms 


0.1 ms 


Idrag 


9 m 


1 cm 


Uflu = 


V Pair 


Pair , 

gd. 



(18) 



Experimentally, the rescaled threshold velocity Qi u is found 
to be around 0.1 (see table [l]). 

Above Ufi u grains spontaneously start rolling at the 
surface of the sand bed p4fl . During their motion, some 
grains can take off the bed, due to the bumps beneath 
them or to the aerodynamic lift force. As soon as they 
have left the bed, the flow around them become again 
symmetric (figure |lj) so that they start being accelerated 
downward by gravity. They collide the bed, rebound and 
eject other grains. The later are accelerated by the wind, 
splash on the sa nd bed and so on, until saturation be 
reached (see part [O). 



3.2.2 Impact threshold 

Suppose that there are already some grains in saltation. 
If the wind velocity decreases, the grains velocity also de- 
creases. So does the sand flux. The minimum wind velocity 
which can sustain the sand transport is that for which the 
grains impact velo city is no more sufficient to eject other 
grains (see section 



3.1.3 ). This happens for Vi mp — a^/gd 
and thus for a wind velocity u* of the order of a few times 
yfgd: 

Uimp — Cirap\/~g~d. (1^) 

Experimentally, the rescaled threshold velocity C,i mp is found 
to be around 3.5 (see table [l]). 



It is somewhat striking that grains can be dislodged by 
the wind since they fall down again after having taken off. 
This is due to the asymmetry of the flow around a static 
grain at the surface of the sand bed, as demonstrated by 
figure |lj. Since r is by definition the horizontal force per 
unit area, the drag force acting on a grain at the sur- 
face of the sand bed scales as Tod 2 . A grain will escape 
from trapping only if this force is larger than gravity. As 
a consequence, there exists a threshold shear velocity u fi u 
below which a sand bed at rest cannot be eroded by the 



3.2.3 Entrainment by the fluid and by impacts 

In the case of aeolian sand transport, we see on table |^ 
that the 'impact' threshold velocity is smaller than the 
'fluid' threshold velocity. This means that the transition 
from the sand bed at rest to the saturated sand trans- 
port presents a hysteresis. Starting from a light wind (it, 
smaller than both Ui mp and Ufi u ), the sand flux is null. 
Then, if the wind velocity increases above the fluid thresh- 
old Ufiu, the grains start being entrained directly by the 
wind and a saturated sand transport establishes. Now, if 
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the wind velocity decreases, the sand flux vanishes at Ui mp 
only. This is characteristic of a subcritical transition. The 
consequences of this hysteretic behaviour has not be in- 
vestigated, so far. 

^From this point of view, the sand transport behaves 
very differently in liquids. It can be seen on table [2] that 
the fluid threshold velocity in water is much smaller than 
the impact threshold. This means that the hysteresis dis- 
appears in water - as in liquids in general. It also indicates 
that the sand transport under water is dominated by di- 
rect fluid entrainment, due to the favourable density ratio. 
By comparison, the splash process is inefficient and there 
is, so to say, no saltation - and therefore no reptation - 
under water. 

In order to distinguish the direct entrainment of grains 
by the fluid from saltation and reptation, we propose the 
name 'tractation' for this motion. This word is derived 
from the latin verb trahere which means to drag Q 'Trac- 
tons' are thus the predominant species under water. Trac- 
tation is very similar to avalanches except that the driving 
mechanism is the fluid drag instead of the gravity. Then, 
the same ki nd o f description could in principle be adapted 
(see section 3.5). 



1 00 



3.3 Saturated flux and saturation length 

3.3.1 Negative feedback on the wind 



Z,, 



ID 



1 



0.1 



0.01- 



m 



2 



u(z) 



Fig. 15. Velocity profile modified by the sand transport, as ob- 
tained theoretically by Owen and Rau pach p3[ and in nu- 
merical simulations by Anderson and Haff 3l||3^ ~It is a piece- 
wise logarithmic profile. Inside the saltation layer (z/Z sa i < 1), 
the shear velocity is decreased but the roughness rd is that of 
the bed. Outside (z > Z 3a i < 1), the shear velocity is u* but 
the apparent roughness zo has increased, as if the soil was 
higher. 



The only limit to the erosion of a sand patch is the satu- 
ration of the sand flux. It occurs because of the negative 
feedback of the transported sand on the wind strength. 
The point is that the same momentum flux is used to 
maintain the turbulent boundary layer and to speed up 
the entrained grains. For a given shear velocity, there is 
thus a maximum flux of sand, called the saturated flux 
Qsati which can be transported by the wind. 

Since they are accelerated by the wind, the grains, at 
a height z, have a smaller velocity Uf when they go up 
than when they come down again (uj.). The volume of 
sand which collides a unit area of the sand bed per unit 
time is, by definition, the flux ip. As shown on figure [l3], 
it is related to the horizontal flux q integrated along the 
vertical and to the hop length I: <j) = q/l. The momen- 
tum transferred to this volume of sand is by definition the 
sand-borne shear stress T san d and is equal to the mass flux 
psand 4> times the velocity difference u-\ — uy. 

Tsand = Psand 4> (l*f ~ lij (20) 

The remaining part of the total shear stress, the air-borne 
shear stress r air = T —T san d, accelerates the flow itself. As- 
suming that the turbulent boundary layer is still at equi- 
librium, the standard turbulent closure (|^) leads a modi- 
fied velocity profile given by: 

1 Saxa ingentia fluctus trahunt (De bello Jugurthino, Sal- 
lustius), i.e. huge stones are dragged by the flow. 



with the same boundary condition than previously, u = 
at the height z — rd. 

Above the saltation layer, there is no grain and thus 
no grain-borne shear stress: r a i r = tq. An undisturbed 
turbulent boundary layer is thus recovered, but with an 
increased apparent roughness zq (figure |l^). Everything 
looks as if the soil has been risen at a height zq — rd. Just 
above the soil, inside the reptation layer, the air-borne 
shear stress r a ir is strongly reduced, and is much smaller 
than tq. As a consequence, the roughness is similar (u 
vanishes at rd) but the apparent shear velocity is smaller 
than (figure p"5| ). Experimentally, the velocity profile 
is found to be logarithmic, with an increased apparent 
roughness Zq, but the lowest part of the profile is too thin 
to be measured [Q,|5[|3^,^,^2). Therefore, the numerical 
and theoretical findings for the velocity profile inside the 
saltation layer should be checked experimentally. 



3.3.2 Equilibrium transport 

In the previous sections, only the consensual properties of 
sand transport have been reviewed. We now come - in this 
subsection - to the controversial points. 

It is clear that the hop height depends on the launch 
velocity w and scales as w 2 jg. The hop length depends 
on w and on the mean horizontal velocity u, and scales 
as wu/g. Now, what are the mean velocity components 
u and w when the equilibrium is achieved? Nalpanis et 
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al. p5[ have found that the mean vertical launch velocity 
w is about 2m* and that both h and I scale on v%/g. On 
the other hand, Rasmussen and Iversen |^,^9|,^2) have re- 
ported measurements of the horizontal grain speed show- 
ing almost no dependence of u on the shear velocity it*. 
The scaling found by Nalpanis & al. |25| is that of the 
saltons while the scaling found by Rasmussen and Iversen 
|}8|,[39[^] corresponds to the reptons. A third scaling has 
been obtained by Jensen and Sorensen @,^7) who have 
used a model to extract the same quantities from experi- 
mental measurements of the vertical variation of the trans- 
port rate. 

Another - indirect - measurement has also been con- 
ducted by Rasmussen and Iversen |38|,[39,42 : they have 



Qsat 



= 



if 



it* < u 



imp 



(24) 



studied systematically the dependence of the apparent 
roughness zq with the wind velocity it* and the grain size 
d. From figure [TJ|, it can be inferred that z is related to 
the height of bounce Z |l3) and to the mean horizontal 
velocity u (~ u(Z)) by zq — Z exp(— u(Z)/u*). If w and 
u scale as it*, zq should scale as it*/<7- If w and it are inde- 
pendent of u* , zq should increase with it* and saturate at a 
value independent of u* . The experimental measurements 
can be rescaled to give approximately: 



-j- + (1 - Qirnpjd 



(22) 



The apparent roughness zq tends to the soil roughness rd 
at the impact threshold and above, increases and scales 
asymptotically as u\j g. According to the previous analy- 
sis, this means that the feedback of the sand transport on 
the wind is dominated by high energy grains. 

As a conclusion, the question of the trajectory mean 
properties is not completely solved, perhaps because it is 
ill posed. It is clear that there is a distribution of charac- 
teristics ranging from the low energy reptons to the high 
energy saltons. Then, the scaling laws obviously depend 
on the way quantities are averaged. For instance, the re- 
sult will be different if an average is weighted by the ver- 
tical flux, by the grain density, by the horizontal flux or 
even by the energy flux. Also the result will depend on the 
height of the lowest trajectory taken into account in the 
statistics, for instance when a camera is used. 



3.3.3 Saturated flux 

The saturation of the sand flux can be understood without 
entering in the details of the mechanisms. The first idea, 
proposed by Bagnold [0 , was that equilibrium is reached 
when the sand-borne shear stress has taken a given part 



of the overall shear stress: T san d cx To 
(f) by q/l in equation (po|), we get: 



Pa 



q sa t oc 



Pa 



I 



T,d Uf 



Replacing 



(23) 



The hop length I scales as uw/g. The velocity difference 
between the rise and the descent — it| should be a frac- 
tion of the grain horizontal velocity it. This gives 



Owen p8[ have introduced a refined argument: the sat- 
uration is reached when the wind shear velocity inside 
the saltation layer has decreased to its threshold value. In 
Owen's article, this threshold was it/; u , meaning that the 
erosion is due - and limited - by the direct aerodynamic 
entrainment. It seems more reasonable to use the impact 
threshold u imp . Since T sand = r - r air = p air u\ - ( imp gd, 
this gives: 



Pair 

Qsat OC W 

Psand 



imp 1 



(25) 



More detailed formulas have been proposed for this rela- 
tion |3^.^,^5) which, as this one, essentially smooth the 
Bagnolds relation around the threshold velocity Ui mp . 

Both expressions ( p4| ) and (25) depend on the typical 
launch velocity w which, as seen in the previous section, 
remains problematic. Nevertheless, we can conclude on the 
scaling of the saturated flux. First, the flux q sat associated 
to high energy saltons (w cx u*) increases more rapidly 
than the flux associated to low energy reptons (w cx \/gd) ■ 
Second, to maintain a possible saturated reptation flux, a 
saltation flux is needed, given by equation (17). For large 
shear velocity, the non dimensional factor l sa i/(lrep N e j e ) 
relating q rep to q sa i scales as u^/^/gd. So, to maintain 
a flux saturated by reptons, a flux of saltons scaling as 
P air I P sandal g is needed. In any case, this is the scaling 
expected at large shear velocity even if reptons can be pre- 
dominant just above the threshold. The last - and the best 
- argument is that aw cx u* is the scaling measured ex- 
perimentally @,^5j,Mj3^], with a prefactor of order unity. 

The expression (|25|) predicts that the saturated flux 
should slightly decrease for an increasing grain diameter 
d, due to the threshold effect. From experimental measure- 
ments, Bagnolds [jlj has reported an increase of q sa t with 
d: for grains ranging from 100 pm to 1 mm, the prefactor 
in equation (|24|) is found to be y/d/D with a new length 
D ~ 150 /xm. This is very striking since D should neither 
depend on u* nor on d. Rasmussen and Iversen ]38|, p9| , |42| | 
have also an extra length in the scaling of the apparent 
roughness zq. Introducing the same length D ~ 150 /im, 
the fit of their data by 



Z ~ ry/D/d (j + 0-- <irnp)<ty 



(26) 



qsat oc 



if u* > u 



imp 



is much better than with expression (p2[). And again, D 
should neither depend on d nor on u*. This new parame- 
ter D have to be related to the fluid entrainment (with a 
Reynolds number dependence) or to the collision process 
(with a dependence on the elastic properties of the mate- 
rial). This implies the existence of new mechanisms, not 
investigated so far. 



3.3.4 Space and time lag 

When the wind arrives at the edge of a sand sheet, it dis- 
lodges and carries away sand grains. The later fall down 
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again and eject other grains so that the transport rate 
increases. In turn, these grains in saltation reacts on the 
turbulent boundary layer to decrease the velocity. There 
is thus a lag between the edge of the sand sheet and the 
achievement of the sand flux saturation. This was first 
reported by Bagnolds H who observed experimentally 
that the transport rate of sand nonetheless increases to- 
wards its saturated value but overshoots and exhibits a 
damped oscillation in space. Other experimental papers 
report space lags but none have recovered the oscillation. 
The existence of a lag was recovered in the numerical sim- 
ulation by Anderson & Haff ]3l] , ^| . Several authors have 
observed that the space lag before saturation - the satu- 
ration length - is almost independent of the shear velocity 
u*, contrarily to the saltation length. Note that in cellular 
automaton models |l4|,[l5|,|l7]], there is also a lag but due 
to transport non locality i.e. to the fact that grains moves 
by jumps of one saltation length l sa i. 
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Fig. 16. Sketch of the promotion of a repton to saltation. 
Once ejected the grain reaches the wind velocity after a typical 
length Idrag- The equilibrium sand transport is achieved when 
the wind speed has so decreased that the number of reptons 
promoted to saltation is just equal to the number of saltons 
absorbed by the sand bed at their impact. 



The physical mechanisms which lead to the sand trans- 
port equilibrium have not to be specified to derive the 
saturated flux. It is the interest of Bagnolds and Owen 
arguments but also their limits. But to understand the 
origin of the saturation length, these mechanisms have to 
be made explicit. Owen pq | have proposed a first idea: 
the flux saturates when the wind does not directly en- 
train grains any more i.e. when the wind shear velocity 
inside the saltation layer has decreased to its threshold 
value Ufi u . It is probably the case for the sand transport 
in water but not in air, in which the direct aerodynamic 
entrainment is much less efficient than the entrainment by 
impacts. 

Anderson & Haff 



32 1 characterise the production 



of new grains in saltation by the mean replacement rate 
77 which is the mean number of grains dislodged by an 
impacting grain (except itself, which rebounds). The flux 
saturates when the wind and thus the grains velocity be- 
comes so small that 77 vanishes. The erosion rate which 
is the difference between the flux of grains taking off the 
bed and the flux of grains impacting on it, is given by the 
spatial derivative of the sand flux d x q sa i- It is equal to 
the incident flux of grain <j) sa i = q S ai/hai times the mean 



a ([sal 
OxQsal = r\- — 



(27) 



The simplest possibility is to imagine, as proposed by 
S0rensen ]29| , that the ejected grains are of the same 
species than the impacting ones. r\ is then simply the num- 
ber of ejected grains N e j e described by the equation (pT[). 
The saturation is reached when N e j e vanishes. This ap- 
pends when all these saltons have a velocity of the order 
of a^/gd - independent of it* - i.e. when they have be- 
come reptons (figure |l5|). Because each grain gives 1+N e j e 
grains after a collision, the flux first increases exponen- 
tially with a typical lengthscale l sa t = l e je equal to the 
saltation length l sa i divided by the typical replacement 
rate r\: 



ad 



gd 



(28) 



This equilibrium situation with a uniform reptation 
layer is far too simplified. If a grain jumps just above this 
layer, it is progressively accelerated and becomes a salton, 
as seen on figure |l6|: the uniform reptation layer is un- 
stable. Second, low energy grains have a large probability 
1 — Preb to be absorbed by the sand bed when they collide 
it (eq. |lj). At equilibrium, this should be balanced by a 
production of reptons. This suggests a slightly different 
picture in which saltons and reptons coexist pacifically. 
The saltons produce reptons when they collide the sand 
bed. reptons are promoted to the rank of saltons, once 
accelerated by the wind (figure nq). This is probably the 
situation reached in Anderson & Haff numerical simula- 
tion ]3l| , p2[ . The production of reptons involves the typical 
length scale l e j e defined above. The acceleration of reptons 
to the velocity of saltons takes a length Idrag scaling on 
the grain size and on the sand to fluid density ratio: 



, ~ Psand j 

''drag — s " 
Pair 



(29) 



As time - and space - goes by, the reptation and salta- 
tion flux increase and the wind strength decreases. Equi- 
librium is reached when the number of reptons promoted 
to saltation just balances the small absorption of grains in 
collisions. Since the process involves two species and thus 
two lengthscales l e j e and Idrag-, the saturation length l sat 
should be the largest of two. It should also be larger than 
the saltation length l sat . For typical values of the shear 
velocity, the promotion of reptons to saltons is the limit- 
ing mechanism so that l sat is given by the inertial length 

Idrag- 

It is interesting to note that there has been only one 
previous derivation of the saturation length |45| which uses 
a completely different argument. Owen's criterion, pre- 
dicts that the air-borne shear stress should have decreased 
to its threshold value at equilibrium. Simultaneously, the 
mean replacement rate rj should vanish. On this basis, 
Sauermann et al. ]45| ] proposed the empirical scaling law 
r\ cc (u*/ufi u ) 2 — 1, which, once reinjected in equation (p7[), 
leads to a saturation length l sat ~ £ p sand I Pair d, as found 
here. 
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Except the logarithmic dependence hidden in £ the 
saturation length is independent of the shear velocity it*. 
It scales as the grain size weighted by the sand to air 
density ratio and it is of the order of 10 m (table |l|), as 
found experimentally fjl|. 



3.4 Coupling between wind and shape 

3.4.1 Erosion rate, minimum size and barchan speed 

In the previous section, we have seen how a sand bed is 
eroded by the wind. Now, where does the erosion takes 
place on a barchan and at which rate? To answer this 
question, let us introduce the local height h(x, y, t) of the 
dune and the volumic sand flux q(x,?/,i). The conserva- 
tion of matter reads: 



<9 t /i + V.q = 



(30) 



If locally the sand flux is larger than the saturated flux 
then the sand flux decreases spatially and sand is deposed. 
If on the contrary the local sand flux is smaller than the 
saturated flux, the sand flux increases spatially and the 
sand bed is eroded. 




Fig. 17. A sand patch (dome) on Tarfaya beach (southern 
Morocco) upwind the barchan field. The precursors to barchans 
directly appear with a length and a width of the order of 10 m. 



What appends if a dune is smaller than the satura- 
tion length I sat? Because it is not saturated, the flux will 
increase continuously over the whole length of the dune. 
The dune will thus be eroded everywhere and will disap- 
pear. This explains very simply the experimental observa- 
tion reported here (figures g and 0) that a small sandpile 
blown by the wind disappears. The order of magnitude 
of l sa t is also the length (and the width) with which the 
sand patches which will give barchans nucleate (figure |i~7| ) . 
The introduction of the saturation length is thus very im- 
portant because it explains the existence of the minimum 
length of the dune. We will see in the subsection devoted 



to the wind around the dune that no other lengthscale 
appears in the problem so that it could be the relevant 
lengthscale in the problem. 

To get a better view of the way a barchan is eroded, we 
can assume that it propagates with a constant shape and 
speed. Then h and q depend only on the variables x — ct 
and y so that the erosion rate dth is equal to —cd x h. It 
means that the dune is eroded in the places where the 
slope along the wind direction is positive and the sand 
is deposed when this slope is on the contrary positive. 
When the brink and the crest coincide (for instance the 
large dune on figure [n]) , the horns and the slip face are 
the only places where the sand accumulates. For a dune 
presenting a small slip face (as the small dune on figure [n]) 
there is a third region of accretion around the brink. 

Let us consider a two-dimensional dune, invariant along 
the transverse direction y direction. Under the hypothesis 
of constant shape and speed, the conservation of matter 
( |30"|) becomes d x (q x — ch) = which immediately inte- 
grates into: 

q(x) =q + ch(x) (31) 

where go is the sand supply i.e. the sand flux on the firm 
soil behind the dune. If the dune if sufficiently long, the 
sand flux q(x) is saturated at the crest of the dune. Then 
the propagation speed immediately derives from the equa- 
tion @: 

qsat — ?0 



H 



(32) 



Now, the three dimensional problem reduces to two di- 
mensions if the conservation of matter (| 30 [ ) i s integrated 
along the transverse direction y. Equation(|3l|) thus holds 
if q and h are replaced by their average across the dune 
width. We make the reasonable assumption that the av- 
erage height along a cross section of the dune scales on 
the overall height H . Then, the barchan speed c scales as 
{qsat — <lo)/H provided that the flux be saturated at the 
crest. This 1/H scaling of the speed was initially proposed 
by Bagnolds For typical values of the shear velocity 
(table [l]), q sat is of the order of a few hundred m 2 /year 
which is th e or der of magnitude found for Q on the field 
(see section 2J5 and in particular equation |lj)). Regarding 
the sand supply dependence, there has been no estimate of 
<7o so far. Field measurements show that Bagnolds scaling 
for the speed is imperfectly verified, a better approxima- 
tion being c ~ Q/(Hq + H). The existence of the cut-off 
scale Hq could come from the fact that the flux at the 
crest is not totally saturated but depends on the ratio of 
the dune length L to the saturation length l sat . 



3.4.2 The wind shape relationship 

By the erosion/ accretion process, the wind modifies the 
topography. But in turn the topography modifies the wind. 
The wind velocity has been measured on the field around 
barchans of heights ranging from 2.5 m to 34 m (fig- 
ure |l^). It turns out that the it is almost independent 
of this height. Even if the vertical variation of the velocity 
is not perfectly logarithmic |47|] , the shear velocity can be 



B. Andreotti, P. Claudin and S. Douady: Selection of barchan shapes and velocities. Part 1 



15 



1.4- 
1 .2- 

i ,0- 
as - 
pi 

0.4- 
0.2 
(1.0 



u*(x)/U 



a) 



h 1 
I 



tin 



if 



h(x)/H 



cS 



0.6 0.4 0.2 0.0 0:2 0,4 0.6 0.8 1.0 

Fig. 18. The white and black circles are measurements of the 
central profile of a barchan (H = 2.5 m and L = 36 m) in the 
Negrita beach (south Morocco). The position x is rescaled by 
the dune length L and the profile h(x) by the dune height 
H. The shear velocity u, on this dune is shown by white 
squares, rescaled by the shear velocity U, far from the dune. 
Three others shear velocity profiles are shown, corresponding 
to barchans of various heights: the black squares (H — 34 m 
and L = 200 m) have been measured by Sauermann et al. [jlfij 
in Jericoacoara (Brazil), the white triangles (H — 10 m and 
L = 100 m) by Wiggs |2fj] in the Sultanate of Oman, and the 
black diamonds (H — 5.5 m and L — 60 m) by Howard et al. 



computed at different place on the dune. At the edge of 
the dune, u* is observed to be slightly smaller than its 
value f/* far from the dune (u*/Z7* ~ 0.9). It then in- 
creases along the dune back: at the inflexion point, u*/C7* 
is around 1.2 and it reaches 1.3 to 1.4 at the top of the 
dune. 

The main theoretical problem is thus to determine the 
flow around a given sand bump and in particular the lee 
recirculations. The scaling laws of the wind velocity field 
are somehow simpler to determine than its spatial vari- 
ations. For a fully turbulent flow, the Reynolds number 
UifH/v based on the dune size, is typically of the order of 
10 6 , and therefore viscous effects are negligible in front of 
inertial effects. This means that the whole velocity field 
is proportional to the shear velocity far from the dune, 
noted [/*. 

Let us first consider the - theoretical - case of a verti- 
cally uniform turbulent flow. An homogeneous turbulent 
flow does not have any proper lengthscale. The whole flow 
is thus scale invariant: two dunes of the same shape but of 
different sizes should be surrounded by the same velocity 
field (in rescaled coordinates). To say it in crude way, if 
the dune size is multiplied by two, the velocity at the top 
remains the same (figure |l^) . 




b) 




velocity 
maximum 



m " 



bump 
maximum 




Fig. 19. A uniform ideal turbulent flow does not have any 
lengthscale, meaning that the velocity at the top of a bump 
(a) is the same if the bump is twice smaller (b). (c) Due to 
a pressure effect, the velocity at the top is larger if the bump 
length is divided by two but not the height. Even if the bump 
presents a symmetry downwind/upwind (d), the velocity field 
does not present the same symmetry: the velocity is larger 
before the crest than after so that the velocity maximum is 
displaced upwind with respect to the bump maximum. 



The main effect is related to changes of pressure [ i8[ 
[l9|| , which is a non local function of the shape. In other 
words, the pressure and thus the velocity at some point 
of the dune is a function of the whole dune shape. Denot- 
ing by h(x) the local height of a two-dimensional dune (as 
that shown on figure fl9| ) , we can conclude from the previ- 
ous considerations that the shear velocity depends only 
on the slope d x h, which is the only dimensionless field de- 
scribing the dune shape. But it is a non local function of 
d x h. Consider the dune shown in figure (a ) and shrink 
it horizontally by a factor of two (figurc[19| c), the wind 
will be larger. This means that the velocity is sensitive to 
the aspect ratio of the bump, to a kind curvature rescaled 
by the size of the dune. This pressure effect is a linear 
effect and exists even in the limit of small bumps. 

If the direction of the wind is a symmetry axis for the 
dune, the wind is symmetrical as well. But if the dune ad- 
mits a symmetry axis perpendicular to the wind (figure |l^ 
d), this symmetry upstream/downstream is broken by the 
velocity field, as in the classical case of the flow around 
a sphere. This corresponds to an irreversible feature of 
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the flow: the streamlines are different when the flow is re- 
versed. This asymmetry upwind/downwind is a non- linear 
effect of the streamlines curvature on turbulence E{| . It is 
an inertial effect which, as previously, is scale invariant. 

If the backward face is steep, the flow can no longer 
follow the form of the dune. The boundary layer separates 
from the bed and reattaches down wind, enclosing a sepa- 
ration bubble. Inside this bubble, the wind flows back to- 
ward the obstacle so that the velocity profile is no longer 
logarithmic. Note that even the exact conditions under 
which flow separates are not clear yet. The best idea to 
model flow separation has been proposed by Zeman and 
Jensen |5Ct| . They suggested that the flow outside the re- 
circulation bubble was precisely that would have been ob- 
served if the streamlines separating the recirculating flow 
from the turbulent boundary layer were solid. Again, this 
idea of a dune envelope seams correct to the first order 
but should be corrected by two effects. First, it is some- 
how arbitrary to prolong the dune by a surface of same 
roughness. Second, the separation surface is well defined 
on the average but fluctuates in time while the actual sand 
bed does not. 

To summarise, the velocity field is scale invariant in 
the ideal case of a vertically uniform turbulent flow. The 
velocity increases in regions where the curvature is nega- 
tive. In the case of a symmetric bump, the velocity field 
around this bump does not present the same symmetry: 
the velocity is larger before the crest than after so that 
the velocity maximum is displaced upwind with respect 
to the bump maximum (figure [n] d) . This means that the 
velocity is also sensitive to the local slope. 

All these properties are present in the expression given 
by Jackson and Hunt [^l],[5^,^3| and simplified by Kroy 
et al. jl9) for the flow around a smooth flat hill. In two 
dimensions, it reads: 



ui(x) = Ut[l + 



A [ — d x h(x 
J tts 



s) + Bd x h(x) (33) 



where A and B are almost constant (see below). Jackson 
and Hunt |5l]]Iiave predicted the values of A and B in the 
limit of a vanishing as pec t ratio H/L - at least smaller 
than 0.05 but equation 33 can also be used for dunes, in- 
troducing effective coefficients A and B. The whole field 
it*(x) is proportional to £/» which is the velocity above 
a flat bed (d x h = 0). It does not depend on the overall 
bump size since it is a - non local - function of the slope, 
only. It is modulated by the local slope (the B term) and 
by a 'dimensionless curvature' (the A term) which takes 
the form of a convolution of the slope by the kernel 1/x. 
Equation (^), used on a dune prolonged by the separation 
bubble modelled empirically, is the best known analytical 
model of wind above a barchan. It has for instance been 
used by Weng et al. Q who computed the erosion rate 
of a barchan and by Kroy et al. |19| who integrated nu- 
merically a complete model of dunes. 

In the case of a turbulent boundary layer, there is ac- 
tually a characteristic lengthscale, the soil roughness zo, 
which breaks the scale invariance. The previous effects re- 
mains valid but the fact that the velocity profile depends 



logarithmically on height leads to logarithmic corrections 
in the dune length to roughness ratio |^1,48,5^,|53|,^9|. 
For instance, in the original Jackson and Hunt model, the 
asymmetry upwind/downwind was directly related to the 
soil roughness (B oc l/ln(L/zo))- We see on this exam- 
ple that the effect of the roughness z is not negligible: 
ln(L/zo) is of the order of 10 for dunes. But it also shows 
that it is almost independent of the dune size: a factor 
10 on L leads to a variation of 20% of ln(L/zo). As a 
conclusion, the wind scale invariance is a very good ap- 
proximation, even with roughness effects. This analysis 
is confirmed by the wind measurements presented on fig- 
ure HI. 



3.5 Avalanches 

The last important phenomenon takes place at the slip 
face. If the sand flux is not saturated upwind the dune, 
the back of the dune is eroded and the sand flux increases. 
This sand is deposed soon after the brink, on the slip face, 
and forms a kind of snowdrift. As shown on figure |20j when 
the slope becomes locally larger than the static friction co- 
efficient an avalanche spontaneously nucleates which 
propagates downward the slip face. It is a dense flow in 
which grains always remain in contact and which is lim- 
ited to a thin layer at the surface of the slip face. As in 
solid friction j54|, this flow stops roughly when the slope 
as decreased below the dynamical friction coefficient /i^. 
The modelling of avalanches have recently received pretty 
high attention from the physicists (5^,^,^|. In particu- 
lar, we now have accurate descriptions of granular surface 
flows by Saint- Venant equations |57],|58|,^9| governing the 
evolution of the free surface, the flowing height and the 
mean velocity. 

Fortunately, we do not need such refined models to un- 
derstand avalanches in the case of dunes. Note first that 
the avalanche duration (a few seconds to one minute) is 
always smaller than the time separating avalanches (a few 
minutes to a few days) which is much smaller than the 
turnover time. Note also that the details of the slip face do 
not react on the dune back since it is inside the separation 
bubble. Each individual avalanche propagates downward 
the steepest slope and stops when the slope has decreased 
below fid- Thus, avalanches may be considered as an in- 
stantaneous slope relaxation process which displaces sand 
along the steepest slope, when the later is larger than jid- 



4 Conclusion: open problems 
4.1 Field observations 

In the first part of this article, we have presented the 
field observations about barchan dunes. Let us give here 
a short summary. A barchan is a crescentic dune propa- 
gating downwind on a firm soil. When the direction of the 
wind is almost constant, these dunes can maintain a nearly 
constant shape and size for very long times. The barchan 
mean velocity c scales approximately with the inverse of 
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Fig. 20. When an avalanche spontaneously nucleates, usually in the middle of the deposition zone, it both propagates downhill 
and uphill and stops when the local slope has decreased below the dynamical angle. 



its height H: c ~ Q/(Hq + H). It turns out that the flux 
Q (typically 100 m 2 /s) and the height cut-off Ho (of the 
order of 1 m) depend on time (probably through meteo- 
rological variations) and on the dune field. The barchan 
height H, its length L, its width W and the length of 



its horns L, u 



are, on the average, related one to the 



others by linear relationships which depend on the dune 
field but apparently not on the wind variations. The typi- 
cal proportions of large dunes are 9 for the ratio width to 
height, 6 for the ratio length to height and 9 for the ratio 
horn length to height. However, H, L and W are not pro- 
portional. It means that barchans are not scale invariant. 
The existence of a characteristic size is confirmed by the 
fact that no barchan lower than 1 to are observed. 

The first thing to note is the strong dispersion of field 
measurements: typically 50 dunes have to be measured to 
establish correctly one relation. As a consequence, there is 
few reliable data despite the numerous studies. Moreover, 
the barchans morphology and speed obviously depend on 
local parameters which have not been identified so far: 
no systematic study has been made on the influence of 
the wind speed, its fluctuations in direction, the nature of 
the soil, the sand supply, the vicinity of other dunes, etc. 
One of the goal, in the future, will be the establishment 
of universal relationships integrating the dependence on 
local parameters. The existence of a minimum size of dune 
also rises questions. What determines this size? If a small 
conical sandpile quickly disappears when eroded by the 
wind, how can barchan form? 

These questions concern barchans in their individual 
behaviour. Another class of problems concerns the global 
dynamics of dune fields (ergs 60 l |6l|). What determines 



the mean spacing between dunes? What selects the size of 
the dune? What determines the sand supply at the back 
of the dune and the leak by the horns? Is sand transport 
more efficient, on the average, in a barchan field when 
compared to saltation over the desert floor? What are the 
precise conditions under which linear transverse dunes ap- 
pear instead of a barchan field? These questions have re- 



ceived pretty low attention. For instance there has been 
one measurement of crescent dune spacing in Namibia by 
Lancaster J62| who have found a linear relation between 
spacing andneight. But it was in a field in which the 
dunes are very close one to the others. Another example is 
provided by Hastenrath (|,^|, who have computed the his- 
togram of dune heights in the Pampa de La Joya (southern 
Peru). In 1964 it exhibits a sharp peak around 3.5 to and 
in 1983 around 1.5 to. These measurements show that the 
dune fields are homogeneous in dune heights but are not 
sufficient to determine the origin of the height selection. 
Concerning the average sand flux, the bulk transport i.e. 
the transport of sand by the dunes have been measured by 
Lettau and Lettau |33| also in the Pampa de La Joya. The 
total flux between 1958 and 1964 was around 1 to 2 /year 
which is a hundred times smaller than the saturated flux 
on a flat soil: bulk transport is apparently much less effi- 
cient than saltation over the flat soil. If the wind is able to 
transport much more sand in saltation, why does the flux 
saturates in the places where barchans form? How and 
where can a dune field form? Moreover they have shown 
that the mean flux increases downwind meaning that the 
soil is eroded (at a rate 200 ^m/year). This article is the 
only tentative of description of sand fluxes at the scale of 
the dune field. 

As a conclusion, further studies will have to focus on 
the detailed characterisation of barchan fields, aiming to 
get sufficient statistics, as for the Pampa de La Joya. 



4.2 Dynamical mechanisms and dune modelling 

In the present state of the art, most of the dynamical 
mechanisms important for barchans formation and prop- 
agation have been identified. The explanation of dune 
propagation is simple: the back of the dune is eroded 
by the wind and the sand transported in the air is de- 
posed at the brink and is redistributed on the slip-face by 
avalanches. The detailed description can be decomposed 
into two parts: the sand transport over a sand bed for a 
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given wind speed on the one hand and the wind speed 
around a given bump on the other hand. 

In this paper, we have proposed a coherent picture of 
sand transport. If the wind strength is sufficiently large, 
grains can be directly entrained by the wind. They roll on 
the soil, take-off, fall back again, rebound, are accelerated 
by the wind, fall back again and so on. Once the transport 
initiated, another mechanism of production takes over the 
direct aerodynamic entrainment. When the grains in salta- 
tion collides the sand bed, they splash up a number of low 
energy reptons which are accelerated by the wind and be- 
come saltons. So, in a first stage, the sand flux increases 
exponentially. To accelerate the grains, the wind has to 
give them some of its momentum. Its strength therefore 
decreases as the sand flux increases. The transport reaches 
equilibrium when the shear velocity has so decreased that 
the number of reptons promoted to saltation just balance 
the number of saltons remaining trapped after collisions 
with the sand bed. This saturation process takes time - 
and space - to establish. The space lag between the edge of 
a sand sheet blown by the wind and the point at which the 
flux saturates is given by the inertial length Idrag defined 
as the distance needed for a repton to become a salton 
i.e. to be accelerated to the wind velocity. This saturation 
length is the only relevant lengthscale of the problem and 
is directly related to the minimum length of barchans. 

There are still a few points concerning sand transport 
which need to be clarified. A first important aim is to re- 
examine experimentally the difference between saltation 
and reptation, in particular the controversial scaling of 
the saltation height and length, the influence of a lateral 
or longitudinal slope, etc. A second problem is the anoma- 
lous scaling of the saturated flux, the apparent roughness 
and probably other quantities, with the grain diameter. 
This suggests the existence of unknown mechanisms which 
would be worth studying. A third interesting study would 
be to reexamine experimentally the saturation process, in 
particular the mechanisms and the scaling of the satura- 
tion length. 

The problem of the turbulent wind flow over a dune is 
different. On the one hand, there exist well known meth- 
ods to solve this problem numerically even though they 
take a lot of computation time. On the other hand, we 
know the basic principles: the velocity field is scale in- 
variant and is proportional to the wind velocity far from 
the obstacle; it increases on a bump and is larger for 
an upwind than for a downwind slope. In middle, very 
few models have been proposed which are both realistic 
and sufficiently simple to be understood. The most use- 
ful is certainly that of Jackson and Hunt The main 
open problem is the description of the recirculation bubble 
which requires either a full 3D simulation or very crude 
empirical assumptions. 

Theoretical and numerical studies of dunes are partic- 
ularly interesting and helpful to understand their dynam- 
ics. The first aim to compute numerically the rate at which 
sand is eroded or deposited on barchan dunes goes back 
to Howard et al. To do so, they used the topography 
of an actual barchan together with laboratory measure- 



ments of the velocity field around a scale model of this 
dune. This work was completed by Weng et al. |fi3f who 
also computed the erosion rate but this time, using Jack- 
son and Hunt approximation. Wipperman and Gross [[l3| 
used a simpler model of flow over hills but computed the 
evolution in time of a sandpilc. They were able to get a 
crescentic dune propagating at a nearly constant velocity. 
However, due to the complexity of their model, this com- 
putation could only be performed on one turnover time. 
More recently, numerical models of the transient of for- 
mation of two dimensional dunes were proposed by Stam 
H and van Dijk et al. @. Finally, Kroy et al. @ have 
proposed a complete model which takes into account all 
the known mechanisms and which also leads to reasonable 
barchan dunes. In part 2, we will simplify this model and 
investigate theoretically the shape and velocity of dunes. 
Much simpler numerical models such as cellular automa- 
tons JlJ,|l5],[l7| have been proposed in the last decade 
which do not lead to realistic isolated dunes but which 
are able to describe patterns of interacting dunes and to 
test the influence of variable winds. The next step is now 
to use numerical simulations or laboratory experiments to 
supply for the lack of control in field studies. They will 
certainly help to find the important parameters for the 
selection of dunes shape, size and velocity but also to look 
at the interaction between dunes and even to test the in- 
fluence of complex wind regimes. 
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